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ABSTRACT 

This work aims to improve the current understanding of the atmospheres of brown dwarfs, especially cold 
ones with spectral type T and Y, whose modeling is a current challenge. Silicate and iron clouds are believed to 
disappear at the photosphere at the L/T transition, but cloudless models fail to reproduce correctly the spectra 
of T dwarfs, advocating for the addition of more physics, e.g. other types of clouds or internal energy transport 
mechanisms. 

We use a one-dimensional (ID) radiative/convective equilibrium code ATMO to investigate this issue. This 
code includes both equilibrium and out-of-equilibrium chemistry and solves consistently the PT structure. 
Included opacity sources are H 2 -H 2 , H 2 -He, H 2 O, CO, CO 2 , CH 4 , NH 3 , K, Na, and TiO, VO if they are present 
in the atmosphere. 

We show that the spectra of Y dwarfs can be accurately reproduced with a cloudless model if vertical mixing 
and NH 3 quenching are taken into account. T dwarf spectra still have some reddening in e.g. 7-// compared to 
cloudless models. This reddening can be reproduced by slightly reducing the temperature gradient in the atmo¬ 
sphere. We propose that this reduction of the stabilizing temperature gradient in these layers, leading to cooler 
structures, is due to the onset of fingering convection, triggered by the destabilizing impact of condensation of 
very thin dust. 

Subject headings: Methods: observational — Methods: numerical — brown dwarfs 


1. INTRODUCTION 

Brown dwarfs with effective temperatures Teff below 
~1000 K (T and Y spectral types) are of great interest to un¬ 
derstand the physics of cool atmospheres and pave the way for 
future studies of cool exoplanets with the James Space Webb 
Telescope (JWST) and the European Extremely Large Tele¬ 
scope (E-ELT). Detailed spectra can now be obtained in the 
near-infrared (NIR) for these objects using e.g. the Hubble 
Space Telescope (HST), or the Gemini Near-Infrared Spec¬ 
trograph. 

The reddening of M/L dwarfs (Teff ^ 2000 K) is thought 
to be the result of the condens ation of silicate dusts and 
the subsequent cloud formation (|Tsuji et al.||1996[ |Chabrier 
|et al.|2000j Allard et al.||2001| l. At the L/T transition (Jeff ~ 
1300 K), this reddening disappears with a sharp transition as 
a function of Teff. Given the hi gh amplitude v ariability ob¬ 
served at the L/T transition (e.g. Radigan|2014 l, this effect is 
interpreted as the clouds breaking up and forming holes at the 
photosphere (e.g.|Marley et al. 12010)1. Yet, the spectral model¬ 
ing of T and Y dwarfs remains a challenge because the IR col¬ 
ors predicted by the cloudless models below ~1000 K appear 
too “blue” compared to observatio ns. This additional redden- 
ing has recently been interpreted by|Morley et al.|(|2012||2014 1 
as the emergence of sulfide and chloride clouds (MnS, Na 2 S, 
and KCl). However, even with this cloudy approach, it ap¬ 
pears difficult to obtain a good model of Y-dwarf spectra (see 
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for example WISEPC J1217b in Leggett et al.|: 

In the present paper, we show under which conditions the 
modeling of T/Y-dwarf spectra is possible with a cloudless 
model. In Sect. we describe our ID numerical code ATMO 
and we apply it to T and Y dwarf atmospheres in Sect. The 
spectral modeling of Y dwarfs is in good agreement w^ ob¬ 
servations when vertical mixing is taken into account. Eor 
T dwarfs, good agreement is obtained when the temperature 
gradient in the atmosphere is reduced. We discuss these re¬ 
sults in Sect. and we suggest that the enhanced cooling in 
the deep layers of T dwarfs arises from the onset of fingering 
convection, resulting from the condensation of some species 
in very thin dust. 

2. METHOD AND CODE DESCRIPTION 

We have developed a ID radiative/convective code ATMO 
solving for the pressure/temperature (PT) structure of an at¬ 
mosphere, assuming a vertical energy balance between the in¬ 
ternal heat flux and the radiative and convective fluxes, and 
vertical hydrostatic equilibrium. The convection is taken into 
account using mixi ng l ength theory (MLT) as decribed in 
Henyey et al.|(|1965|) andjGustafsson et al.|(|20081) with a mix- 
ing length of 1.5 tirnes the local pressure scale height. 

The PT structure is solved on a cartesian plane-parallel 
grid for a given effective temperature and surface grav¬ 
ity log g. The chemistry takes into account ~150 species at 
equilibrium with the condensation of ~30 liquids/solid s as¬ 
suming elemental abundances from|Caffau et al.|(|201 l|l, and 
is solved by minimization of Gibbs free energy! The refer¬ 
ences of the th ermodynamic data for the gas-phase species 
can be found in |Venot et al.| ( 20 T^ except for H 2 whose ther¬ 
modynamic data were htted from the S aumon-Chabrier equa¬ 
tion of state (see|Saumon et al.|1995|l. The data for the con 


densed species can be found in |Robie and Waldbaum (1968 1 
and|Chase|(|1998|). In the present models, we neglect ram-out 
processes and clouds. This approximation is critical to as- 
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sess the abundances of some refractory elements such as Na, 
K, and S because of the different condensation paths that can 
happen depending on the thermal history of the atmosphere 
(e.g. NaAlSi 308 /KAlSi 308 will not form if A1 and Si are re¬ 
moved at higher temperatures and Na and K will condense 
in Na 2 S and KCl). However in the spectral window of in¬ 
terest, the main absorbers such as H 2 O, NH 3 , CH 4 , CO are 
not affected by this hypothesi s and their abundances can even 
be computed analytically (see |Burrows and Sharp 1999|l. The 
code is coupled to the kinetic chemical network of Venot et ah] 
( | 2012 [ ), which includes ~1000 reversible reactions associated 
with ~100 species based on C, H, N, O elements. Vertical 
mixing is taken into account with a free parameter repre- 


senting the vertical eddy diffusion coefficient (see Ackerman 
and Marley|2001 1 . We perform a full chemical network run on 
the atmosphere by integrating the chemistry in time (usually 
up to 10 '® seconds) and we reconverged the PT structure of 
the atmosphere on the fly until both the chemistry and the PT 
structure are in a steady state. The PT structure can be solved 
consistently with the out-of-equilibrium chemistry and can 
have a significant difference if a main absorber is quenched 
below the photosphere by the vertical mixing. We test ed the 
kinetic chem ical network by reproducing the results of |Venot| 


|etar] ( |20T2l l for HD209458b and HD189733b. Our result Tare 


also in good agreement withjHubeny and Burrows |(j2007|), es¬ 
pecially considering the differences in elemental abundances 
and opacity line lists. For our case at 7’efr=1000 K, log g=4.5 
and log /rzz=8.0, we predict a quenched abundance of CO 
at ~2xl0 and of NH 3 at ~7xl0“^, an d the model with 
r..ff-900 K, log g=4.5 and log in|Hubeny and Bur^ 


rows 


(2007 I indicates a CO abundance of ~ 1.5 x 10 (fast2 


chemical timescale|Yung et al.|1988|l and aNH 3 abundance of 
~ 5 X lO-'’. 

The line-by-line radiative transfer has a lready been de¬ 
scribed and tested in |Amundsen et al[] ( |20141 l. The differences 
in radiative flux and heating rate between our line by line code 
and the radiative scheme used in the global circulation model 
of the UK Met-Office, the Unifed Model ( [Amundsen et al.| 
2014|l are generally of the order of a few percent. We use the 
same opacities for H 2 -H 2 , H 2 -He, NH 3 , H 2 O, CO, TiO and 
VO (see Table 1 in [Amundsen et aL]|2014| l. Methane (CH 4 ) 
has bee n updated with the new line li st from the Exomol 
project ([Yurchenko and Tennyson|2014|l and we included the 


CO 2 o pacity using the line list from Tashkun and Perevalov 
( 201 1 || with the line broadening available in the literature 


(broadening coefficients fr om Thibault et al.[ (1992 2000 1 ; 
Sharp and Burrows[([2007[l; [Padmanabhan et al.[([2014|l). We 
have also calculated opacities due to the extremely pressure- 
broade ned sodium and pota ssium doublets using the line pro¬ 
files of [Allard et al.[ ([2007[l. We have used the correlated-A: 
method described in[Amundsen et al.[(|2014|l to speed up the 
convergence and conibined the fe-coe mcients by using the ran - 
dom overlap method for the mixture ( Lacis and Oinas|1991| l. 
Then we compute the surface spectrum from the converged 
PT structure using the line by line radiative transfer. The ra¬ 
diative transfer module also includes Rayleigh scattering by 
H 2 and H e, using the accelerated A-i teration technique de¬ 
scribed in[Bendicho and Bueno[([T995]l. 


3. APPLICATION TO Y AND T DWARFS 
3.1. Y dwarfs and NH3 quenching 
Saumon et al. ( 2006|l have already demonstrated that 


dwarfs based on the observed strength of th 9-11 /rm absorp¬ 
tion feature, ft has also later been confi rmed in the NfR lines 
of NH 3 (e.g. Saumon et al.[[201^ , and Cushing et al. ( [2011 1 
confirmed the signature of vertical mixing in all but one of 
their Y dwarfs by the very absence of more detectable NH 3 
absorption, ft has been expected that absorption by ammonia 
in the Y band should be present, but so far the data have no t 
identified such a signature (e.g. Fig. 9 in Cushing et al.[201 1 1 , 
ft has been also recently proposed by Leggett et al. ( [2015[ l that 
ammonia is depleted in the atmospheres of Y dwarfs. Using 
the Y dwarf WfSEPC J1217b (YO, [Leggett et al.|2014| l as a 
template, we show in the left panels of Fig. f^that indeed a low 
abundance of ammonia is required to reproduce the shape of 
the Y band flux around 1 yum and the correct flux in the H 
band. Indeed, as shown in Fig. when vertical mixing is 
taken into account with a diffusion coefficient, the abundance 
profile of ammonia is quenched at a relatively deep level (50- 
100 bars), resulting in a depletion of ammonia by approxima- 
tively a factor 3, significantly reducing the absorption in the 
peaks in the Y band. We took a value of - 10* cm^s ' 
for the eddy diffu sion coefficient, which is wi thin the high 
values explored by [Hubeny and Burrows[ ( [?007| l; however the 
quenching of the NH 3 /N 2 chemistry is relativ ely insensitive 
to the choice of (see Cushing et al.|20lT| l. When we in¬ 
clude the vertical mixing, the condensation of H 2 O and NH 3 
is ignored for simplicity, this approximation is not an issue be¬ 
cause the difference with and without condensation at equilib¬ 
rium is small at these effective temperatures (Teff is shifted by 
~10 K). Importantly for the color-magnitude diagrams in the 
IR, the flux in the H band at 1 .6 jum is increased and matches 
well the observed spectrum. As shown in Fig.|T] this increase 
in the flux is sufficient to produce the reddening required in 
J - H to match the observations. Applying the same physical 
treatment to the coo ler brown dwarf WISEPA J1541 (Y0.5, 
Cushing et al.|20lT] ), we obtain the same significant improve- 
ment between observed and synthetic spectra, as shown in the 
bottom-left panel of Fig. [Tj Si nce we have included the com¬ 
plete chemical network of[Ven 
predict the quenching of 


impact the flux at 4.5 yum. 


vfenot et al.[( [2012[) in ATMO, we also 
(see Fig.[^, which will 


3.2. Tdwarfs and reduced temperature gradient 

T dwarfs have bluer colors in J-H and J-K than L dwarfs 
which indicates that the silicate and iron clouds should fall 
below the photosphere at the L/T transition. However, their 
spectra still present a reddening compared to cloudless models 


and [Morley et al.[ 
such as sulhde cl 


proposed that another type of clouds 
ouds might be responsible for this residual 


monia is well below its chemical equilibrium abundance in T 


redde ning. We have calculat ed spectra for the T d warfs ROSS 
458C ([Burgasser et al.|2010[l and UGPS 0722-05 ([Lucas et al.| 
2010j l whose spectra ^e shown in rig ht panels of Fig. As 
suggested by [Burgasser et al.[ ( [2010| i, we used a supersolar 
metallicity ([My/H]=0.3) for ROSS 458C based on the mea¬ 
sured metallicity of its two companions ROSS 458 A and B. 
Our cloudless models overestimate the flux in the Y and J 
bands, yielding significantly bluer colors in 7 - HjJ - K com¬ 
pared to observations. Since the photosphere in the Y and 
J bands is deeper than the the photosphere in the H and K 
bands, reducing the temperature gradient in the atmosphere 
should decrease the YjJ flux compared to the HjK flux. We 
have tested this idea by artificially modifying the adiabatic in¬ 
dex y - CpfCv of the atmosphere. The pressure/temperature 
profiles obtained with y-1.2 for ROSS 458C and 7=1.27 for 
UGPS 0722-05, respectively, are displayed in Fig.|^ and the 
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Figure 1. Left panel, top: spectral flux of WISEPC J121756.91+162640.2b (YO) in red compared to a cloudless model with no vertical mixing (green) and 
with vertical mixing - 10^ cm^s“^, blue). Lef t panel, bottom: similar m odels for the sp ectral flux of WISE PA J15415L66225025.2 (Y0.5). The distances 

respectively. Right panel, top: spectral flux of 


are taken from parallax measurements: 10.1^]-^ pc Dupuy and Kraus 2013 and 5.7^® !^ pc i Tinney et al 


2014 


ROSS 458C in red compared to a cloudless model (green) and with a reduced temperature gradient (blue). Ki^t panel, bottom: similar models for the spectral 


flux of UGPS J072227.51-054031.2. The distances are taken from parallax measurements: pc van Leeuwen 2007 and 4.1^^-g pc iLucas et al. 

respectively. All the models assume a radius of 0.1 R©, no adjustment of the total flux is used. 


2010 , 



Abundances 


Figure 2. Abundance profiles of the main species in the modeling of 
WISEPA J1541. The dashed profiles coiTespond to the model at chemi¬ 
cal equilibrium and the solid ones to the model with vertical mixing with 

corresponding spectra in Fig. [T] These models correctly re¬ 
produce the spectra and the colors in 7 - HjJ - K without the 
need to invoke additional cloud opacity. As shown in Fig. 
the amplitude of the temperature reduction required to matdi 
the observed spectra decreases with the effective temperature. 
Therefore, for Y dwarfs we found that there is no need to 



Figure 3. Pressure/temperature profiles of the models of ROSS 458C and 
U(TPS 0722-05. The solid lines con'espond to models with modified adia¬ 
batic index y. The solid circles indicate the minimum pressures for which 
convection is present. Without a modified adiabatic index, the convection is 
overturning, while the modified index, illustrates heat transport by the finger¬ 
ing convection. The change in Kzz has little impact on the PT profile. 

adjust the adiabatic index with respect to its expected value 
(~ 1.45). The relative importance of the effects of vertical 
mixing and the reduction of the temperature gradient can be 
seen on the color-magnitude diagrams in Fig. 
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4. DISCUSSION: ATMOSPHERIC FINGERING CONVECTION 


As shown in Sect. 3.1 a coherent cloudless model taking 
into account vertical mixing in the atmosphere is able to re¬ 
produce the near-infrared spectrum of Y dwarfs. Figure ffl 
shows the color-magnitude diagrams for various cool T and 
Y brown dwarfs in My vs 7 - // and My vs Y - J. We have 
plotted a series of different models at constant log g = 4.5 as 
well as the c loudless and sulfide/H 20 cloud models of|Morle^ 


et al.|(|20141l. We display the models with and without vertical 
mixing and with and without modified adiabatic index to il¬ 
lustrate the relative impacts of these processes. Although the 
cloudy models do yield redder colors in J-H, pr oviding a bet¬ 
ter ma tch to the observed colors (see Fig. 16 in|Morley et al. 


( |2014[ )), they predict also a large reddening in Y - 7 which 
is incompatible with the observations. The problem is even 
worse at larger gravity (log g = 5). In contrast, our cloud¬ 
less models reproduce the correct J-H and Y - J colors for 
Y dwarfs, with out any adjustable parameter. Furthermore, as 
shown in Sect. 3.1 cloudless models consistently taking into 
account departure from chemical equilibrium due to vertical 
mixing in the atmosphere nicely reproduce the near-infrared 
spectrum of Y dwarfs. Our cloudl ess models are shifted in 
J-H with respect to the ones of |Morley et al.| (|2014|l by 
nearly one magnitude at log g=4.5 and we would get similar 
offsets if we compare our models at log g=4.5 to their models 
at log g=5.0. Part of the shift (0.3 mag) can be explained by 
changing in our calculations, (i) the elemental abundances to 
[Fodders (|2003|l, (ii) the line list for methane to the STDS line 
list ( Wenger and Champion||T998) l, (ii) the thermodyna mical 
data of H 2 to the one of the JANAF database as used in jMor-j 
ley et al. ( 2014| l. The origin of the rest of the shift remains 
unknown but is probably due to differences in the opacities 
of the other molecules. It seems more likely, however, that 
such low-temperature low-mass objects have gravities close t o 
log g=4.5 rather than log g=5.0 (IChabrier and Baraffe|2000||. 

For T dwarfs, both our cloudless models and the cloudy 
models reproduce the 7 - H and Y - J colors. In the for¬ 
mer case, as shown in Sect. |3.1| a good match of the spec¬ 
trum requires a lower value of the atmospheric characteristic 
adiabatic index (a good match is also obtained in My versus 
7 - K). This adjustment can thus reflect either the effect of 
clouds or the signature of an energy transport mechanism that 
reduces the temperature gradient in the atmosphere. Given the 
fact that, for T dwarfs, iron and silicate clouds are expected to 
have settled below the phostosphere, as mentioned previously, 
we favor the second suggestion. We suggest that the cooling 
is due to the onset of fingering (thermohaline) convection in 
T dwarf deep atmospheric layers, simil ar to the one trigg ered 
by salt gradients in Earth oceans (e.g. Rahmstorf |2003| l. In 
the limit of “infinitely-thin” dust condensation, we have com- 
puted the mean-molecular-weight gradient in the atmosphere, 
Vj^ assuming that all the condensates remain in the atmo¬ 
sphere in their “mono-molecular” form (e.g., assuming that 
equilibrium chemistry predicts a given abundance of solid 
MgSiOs, we compute the mean molecular weight by keep¬ 
ing this abundance with the molecular weight of MgSiOs). 
This mean-molecular-weight gradient should be fairly repre¬ 
sentative of in the presence of small grains if the grain 
formation process does not si gnificantly change throug hout 
the atmosphere. As shown injRosenblum et al.|(|2011|l (see 
also [Leconte and Chabrier||2012jl, the value of the dimen- 


> Vx = dln(X)/dln(P) 


sionless parameter Rq - (Vy - Vad)/V^ determines the ex¬ 
tent of the overturning convection zone {Rq < 1). In case 
the mean molecular weight increases with height in a sta¬ 
ble atmosphere (V^ < 0, and (Vy - Vad) < 0), as in the 
present context, the possible existence of fingering convec¬ 
tion is determined by the relation 1 < Rq < 1 / t , where 
T = k^Ikt defines the inverse Lewis number, i.e. the ratio of 
the molecular diffusivity to the radiative thermal diffusivity 


kt - \6crT^liJiKp^Cp). This is the opposite situation to oscil¬ 
latory double-diffusive convection, due to a positive molecu¬ 
lar weight gradient and a destabilizing temperat ure gradient, 
suggested to occur in some giant planet interiors ([Leconte and| 
Chabrier[[2012[ [2013[|. Figure [S] shows the atmospheric pro¬ 
files of 17 (tRo) f^OSS 458C and WISEPA J1541, taking 


into account the molecular diffusivity of MgSi 03 (s) and as¬ 
suming that the thermal diffusivity is dominated by the ra¬ 
diative processes {kj varies from ~ 10 “* cm^s“* at 1 kbar to 
~ 3 X 10^ cm^s“* at 1 bar). As seen, fingering convection 
(i.e. 1 /(tRo) > 1 ) is predicted to occur in some parts of these 
atmospheres. The magnitude of the mean-molecular-weight 
gradient is modest; |V^| ~ lO^^^ - lO^^ for ROSS 458C. Inter¬ 
estingly, while the extension of the fingering-convection zone 
is very large for ROSS 458C, it remains relatively small for 
the Y dwarf WISEPA J154L This is fully consistent with the 
fact that we need a lower adiabatic index for T dwarfs whereas 
this modification is not necessary for Y dwarfs. 

The transport mechanism needed to lower the temperature 
gradient has to be efficient to transport energies at the level 
of the internal flux of the object. The possibility that finger¬ 
ing convection could carry out such energy fluxes remains to 
be demonstrated. Also the mixing induced by the convective 
fingers will reduce the destabilizing mean-molecular-weight 
gradient. A consistent theory of fingering convection in our 
model is needed in order to assess the possibility to have a 
self-sustained steady-state energy flux from the convection. 
Nevertheless our results strongly suggest that fingering con¬ 
vection could occur in T dwarf atmospheres due to the desta¬ 
bilizing effect of dust condensation in the stably stratified part 
of the atmosphere, potentially impacting the energy transport 
in these objects compared to both warmer and cooler objects 
and leading to efficient cooling of their deep atmospheric lay¬ 
ers. It can be expected that the increase of the convective 
energy transport will reduce the radiative flux in the corre¬ 
sponding region, hence reducing the temperature gradient in 
the atmosphere, as needed for the reddening in 7 - //. 


5. CONCLUSIONS 

In the present paper, we have shown that 


• Cloudless atmosphere models can reproduce the ob¬ 
served infrared Y dwarf spectral energy distribution, 
provided vertical mixing and out-of-equilibrium chem¬ 
istry are properly taken into account to correctly predict 
the quenching of ammonia at deep levels, modifying its 
abundance profile. 

• Cloudless models for T dwarfs with a reduced tempera¬ 
ture gradient in the atmosphere correctly reproduce ob¬ 
served fluxes and colors. Such a reduced gradient can 
be obtained by a modification of the adiabatic index 
in the atmosphere and could reflect either the effect of 
clouds or of another type of energy transport in the at¬ 
mosphere. If clouds are effectively responsible for the 
reddening, the modification of the adiabatic index is an 
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logg=4.5 



__ logg=4.5 _ 

- sulfide/H20 clouds, Morley et al?t2014) 

- - cloudless, Morley et al. (2014) 

cloudless, logkzz=0, modifiec^- 

- - cloudless, logkzz=0 ^ 



Y-J [MKO] 


Figure 4. Color magnitude diagrams (top: J-mag vs J - H, bottom: Y- 
mag vs F - /) for T/Y brown dwarfs. The Y-dwarf photometi^ is from 


Dupuv and Kraus 

2013L IBeichman 

et al. 

^2014) and the L/T/M from 

Dupuy and Liu 

■p 

jl2); lEaherty et at. 


. Ihe modified adiabatic in- 

dex IS 'y=(L15, 

1XE23) tor 7eff=(l0U0K:8DDK,600K), no modification for 


Yeff <600K. 

easy way to mimic the effect and could be used to better 
constrain cloud models that are currently used. 

• We suggest that fingering convection could be respon¬ 
sible for such a reduced temperature gradient. We 
demonstrate that the condensation of very thin dust un¬ 
der typical T dwarf atmosphere conditions could trig¬ 
ger this mechanism and that the extent of this process 
decreases with decreasing effective temperature, essen¬ 
tially vanishing for Y dwarf atmosphere conditions. 

Fundamental physical mechanisms such as atomic diffu¬ 
sion and hydrodynamical instability might thus take place in 
cold brown dwarf atmospheres and play a major role in their 
spectral evolution. Future observations of cool objects with 
SPHERE, GPI, and JWST combined with comparisons with 
the different models should enable us to distinguish between 
these two effects: presence of clouds or reduced temperature 
gradient in T dwarf atmospheres. Our ability to constrain 
which physics is indeed present in cool brown dwarf atmo¬ 
spheres will bear important consequences for the future un¬ 
derstanding of cool exoplanet spectra. 
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Figure 5. Profile of l/(rRo) for ROSS 458C and WISEPA J1541. Fingering 
(thermohaline) convection is possible where 1/{tRq) > 1. The dotted hori¬ 
zontal line shows the limit between the overturning and fingering convection 
zones for ROSS 458C. 
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